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Abstract. — Fractional Pontryagin's systems emerge in the study of a class of fractional optimal 
control problems, see [2l [3J [9l I17L 1181 124| and references therein, but they are not resolvable 
QQ . in most cases. In this paper, we suggest a numerical approach for these fractional systems. 

Precisely, we construct a variational integrator allowing to preserve at the discrete level their 
intrinsic variational structure. The variational integrator obtained is then called shifted discrete 
fractional Pontryagin's system. 

In [S], we have provided a solved fractional example in a certain sense. It allows us to test 
in this paper the convergence of the variational integrator constructed. Finally, we also provide 
a discrete fractional Noether's theorem giving the existence of an explicit computable discrete 
constant of motion for shifted discrete fractional Pontryagin's systems admitting a discrete sym- 
metry. 
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Introduction 

m 

^^ ■ The fractional calculus, i.e. the mathematical field dealing with the generalization of the 

derivative to any real order, plays an increasing role in many varied domains as economy [12] 
or probability |26L 135] . Fractional derivatives also appear in many fields of Physics (see [22] ) : 
wave mechanic [5], viscoelasticity [6], thermodynamics [23], fluid mechanic in heterogeneous 

S^ | media |20L 137] 138] . etc. Recently a subtopic of the fractional calculus gains importance: it 

concerns the variational principles on functionals involving fractional derivatives. This leads 
to the statement of fractional Euler-Lagrange equations, see PQ, [§l 133] . 

A direct consequence is the emergence of works concerning a particular class of fractional 
optimal control problems, see [21 El El 117], 118] 124] and references therein. Using a Lagrange 
multiplier technique or not, authors obtain with a calculus of variations a necessary condition 
for the existence of an optimal control. This condition is commonly given as the existence of a 
solution of a system of fractional differential equations called fractional Pontryagin's system. 
Hence, the explicit computation of a potential optimal control, from the above necessary 
condition, needs the resolution of a fractional Pontryagin's system which is a main drawback. 
Indeed, solving a fractional differential equation is in general very difficult. Moreover, a 
fractional Pontryagin's system involves left and right fractional derivatives which is an 
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additional obstruction. 

In this paper, we then develop a numerical approach. Let us remind that there exist many 
works concerning the statement of discrete operators approaching the fractional derivatives 
(see |15L 1191 131] ) and then concerning numerical schemes for fractional differential equations 
(see [14|, I27|, I30L 129] ). In particular, one can find studies concerning the discretization of 
fractional Euler-Lagrange equations [4], lllj and fractional Pontryagin's systems [2], SI El H3|, 



Nevertheless, a fractional Pontryagin's system admits an intrinsic variational structure: 
its solutions correspond to the critical points of a cost functional. Moreover, this variational 
structure induces strong constraints on the qualitative behaviour of the solutions and it 
seems then important to preserve it at the discrete level. A variational integrator is a 
numerical scheme preserving the variational structure of a system at the discrete level. We 
refer to Section [2] for more details concerning the construction of a variational integrator 
and let us remind that the variational integrators are well-developed in |21L |28| for classical 
Euler-Lagrange equations and in [11] for fractional ones. In this chapter, we construct a 
variational integrator for fractional Pontryagin's systems and it is called shifted discrete 
fractional Pontryagin's system. 

In [9], we have suggested a deviously way in order to get informations on the solutions of 
a not solvable fractional Pontryagin's system. Indeed, we have stated a fractional Noether's 
theorem giving an explicit constant of motion for fractional Pontryagin's systems admitting 
a symmetry. We refer to [9j for more details and we remind that this result is based on a 
preliminary result proved by Torres and Frederico in [17|, 118] . In this paper, following the 
strategy of the continuous case, we introduce the notion of a discrete symmetry for shifted 
discrete fractional Pontryagin's systems and we finally provide a discrete fractional Noether's 
theorem giving an explicit computable constant of motion. 

The paper is organized as follows. Section [1] is devoted to a reminder on the fractional 
calculus and on the emergence of fractional Pontryagin's systems in the study of a class 
of fractional optimal control problems. In Section [21 after a reminder concerning discrete 
fractional derivatives, we focus on the construction of a variational integrator for fractional 
Pontryagin's systems. We make some numerical tests in Section [3j Especially, let us remind 
that a fractional example is solved in [9] in a certain sense. Consequently, we can test the 
convergence of the variational integrator both in the classical and strict fractional cases. 
Finally, Section 0] is devoted to the statement of a discrete fractional Noether's theorem. 
Technical proofs of Lemmas are provided in Appendix [A! 



1. Reminder about fractional Pontryagin's system 

In this section, we first make a reminder about fractional calculus in Section 11.11 Then, in 
Section 11.21 we remind how fractional Pontryagin's systems emerge from the study of a class 
of fractional optimal control problems. Let us introduce the following notations available in 
the whole paper. Let a < b be two reals, let d, m G N* denote two dimensions and let || • || be 
the euclidean norm of M and R m . 



1.1. Fractional operators of Riemann-Liouville and Caputo. — The fractional 
calculus concerns the extension of the usual notion of derivative from non-negative integer 
orders to any real order. Since 1695, numerous notions of fractional derivatives emerge over 
the year, see |25L I32L I34j . In this paper, we only use the notions of fractional integrals and 
derivatives in the sense of Riemann-Liouville (1847) and Caputo (1967) whose definitions are 
recalled in this section. We refer to I25L I32L 1341 for more details. 



Let g G ^°([o, 6],M ) and a > 0. The left (resp. right) fractional integral in the sense of 
Riemann-Liouville with inferior limit a (resp. superior limit b) of order a of g is defined by: 



t 

Q!-l 



Vt £]o, b], rg{t) := — - / (t - y) a ~ L g(y) dy (1) 

r («) Ja 

respectively: 

Vt G [a, b[, I$g(t) := -L / (y - i)*" 1 ^) dy, (2) 

where T denotes the Euler's Gamma function. For a = 0, let l£g = FLg = g. 

Now, let us consider < a < 1. The left (resp. right) fractional derivative in the sense of 
Riemann-Liouville with inferior limit a (resp. superior limit b) of order a of g is then given 
by: 

Vt G]o,6], Dlg{t) := |(rt)(*) (resp. Vt G [a, b[, D a + g{t) := - 1 ( I^g) (t)) , (3) 
provided that the right side terms are defined. 

In the Riemann-Liouville sense, the strict fractional derivative of a constant is not zero. 
Caputo then suggests the following definition. For < a < 1, the left (resp. right) fractional 
derivative in the sense of Caputo with inferior limit a (resp. superior limit b) of order a of g 
is given by: 

VtG]a,6], c Dlg{t):=Dt{g-g{a)){t) (resp. Vt G [a, b[, c D<lg{t) := D%(g - 5 (6)) (t)) , 

(4) 
provided that the right side terms are defined. Let us note that if g(a) = (resp. g(b) = 0), 
then c D°g = Dig (resp. c D«g = D%g). 

In the classical case a = 1, the fractional derivatives of Riemann-Liouville and Caputo both 
coincide with the classical derivative. Precisely, modulo a (—1) term in the right case, we 
have Di = c Dl = -D\ = - C D\ = d/dt. 

1.2. Reminder about a class of fractional optimal control problems. — From now 
and for all the rest of the paper, we consider < a < 1 and A G R d . Let us denote [a] the 
floor of a. 

In this section, let us remind the following definitions concerning the class of fractional optimal 
control problems studied in [9]: 

— The elements denoted u G ^°([a, 6],M m ) are called controls; 
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- Let / be a ^ 2 function of the form: 

/ : R d xR m x [a, b] — ► R d (5) 

(x,v,t) i — > f{x,v,t). 

It is commonly called the constraint function. We assume that / satisfies the following 
Lipschitz type condition. There exists M > such that: 

V(x 1 ,x 2 ,v,t)£(R d ) 2 xR m x [a,b], \\f(x 1 ,v,t)-f(x 2 ,v,t)\\ < M\\x x - x 2 ||; (/ x lip) 

- For any control u, let q u,a G ^' a '([a, 6],R ) denote the unique global solution of the 
following fractional Cauchy problem: 



,D°_q = f(q,u,t) , ■ 

q(a) = A. { ^ i 

qU,a - g comm only called the state variable associated to u. Its existence and its unique 



ness are provided in [9] from Condition (f x lip) 



— Finally, the fractional optimal control problem studied in [9] is the problem of optimiza- 
tion of the following cost functional: 

C a : tf°([a,&],R m ) — »• R (6) 

M »— »• /" L(q u ' a ,u,t)dt, 

J a 

where L is a Lagrangian, i.e. a ^ 2 application of the form: 

L : M d x R m x [a, b] — ► R (7) 

(x, v,t) i — )• L(x,v,t). 

A control optimizing C a is called optimal control. A necessary condition for a control u to be 
optimal is to be a critical point of C a , i.e. to satisfy: 

VuG^°([a,&],R m ), DC a {u){u) := lim ^ + £U) - C ° {u) = 0. (8) 

£->o e 

In [9] , we then focused on the characterization of the critical points of C a . Firstly, we proved 
with an usual calculus of variations the following Lemma [1] giving explicitly the value of the 
Gateaux derivative of C a : 

Lemma 1. -- Let u, u G ^°([a, 6],R m ). Then, the following equality holds: 

— ( q ^,u,t).q+—(q u > a ,u,t)-udt, (9) 

where q G c rf^ a '([a,b],W i ) is the unique global solution of the following linearised Cauchy 
problem: 

c D«q = ^-(q u > a ,u,t) x q + ^f(q u > a ,u,t) x u 

ox dv (LCP^) 

q(a) = 0. 

This last result not leading to a characterization of the critical points of C a , we then introduced 
the following elements stemming from the Lagrange multiplier technique: 



Let H be the following application 

H : R d x R m x R d x [a, b] 



(10) 



(x,v,w,t) i — ► L(x,v,t) + w ■ f(x,v,t). 

H is commonly called the Hamiltonian associated to the Lagrangian L and the constraint 
function /; 

For any control u, let p u,a G ffl a >([a, b],R d ) denote the unique global solution of the 
following fractional Cauchy problem: 

T 



f)J-f 

;D%p=—(q u ' a ,u,p,t) 



dL t of 

— (q^,u,t)+l-^(q u > a ,u,t) 



(cp^; 



p(b) = 0. 

pU,a - g C ommonly called the adjoint variable associated to u. Its existence and its 
uniqueness are also provided in [9]. Let us note that c D®p u ' a = D"p u ' a since p u ' a (b) = 0. 

Consequently, for any control u, the couple (q u ' a ,p u ' a ) is solution of the following fractional 
Hamiltonian system: 

r)ff 
zD°Lq = -^—(q,u,p,t) 

(HS Q ) 



dw 



D%p=—(q,u,p,t). 

Finally, the introduction of these last elements allowed us to prove the following theorem: 

Theorem 2. - - Let u G ^°([a, 6],lR' m ). Then, u is a critical point of C a if and only if 
(q u ' a ,u,p u,a ) is solution of the following fractional stationary equation: 



9H, , n 

— {q,u,p,t) =0. 



(SE C 



From Theorem [21 we retrieved in [9] the following result leading to the fractional Pontryagin's 
system: 

Corollary 3. C a has a critical point in ^°([a, 6],R m ) if and only if there exists (q,u,p) G 
<af H ([a, b],R d ) x ^°([a, b], R m ) x ^W ([a, 6], K d ) solution of the following fractional Pontryagin's 
system: 

r)ff 
c D"q= —(q,u,p,t) 



flu 

D%p=-^(q,u,p,t) 

8H, 

— (q,u,p,t) = 



(PS C 



[ {q(a),p(b))=(A,0). 
In the affirmative case, u is a critical point of C a and we have (q,p) = (q u ' a ,p u ' a ). 

Let us note that the fractional Pontryagin's system (|PS a p is made up of the fractional Hamil- 
tonian system (|HS Q p . the fractional stationary equation (|SE a p and initial and final conditions. 



In practice, see Examples in [9], we use more Corollary [3] than Theorem [21 Let us remind 
that Corollary [3] was already provided in (2J [31 I17|, I18|, I24| and references therein without 
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Condition (f x lip). However, this result is proved, in each of these papers, using a Lagrange 



multiplier technique requiring the introduction of an augmented functional. In [9] , Condition 



( f x lip ) allowed us to give a complete proof of this result using only classical mathematical 
tools adapted to the fractional case: calculus of variations, Gronwall's Lemma, Cauchy- 
Lipschitz Theorem and stability under perturbations of differential equations. We refer to [9] 
for more details and for a discussion on the subject. 

As we have seen in this section, fractional Pontryagin's systems emerge from the study 
of a class of fractional optimal control problems. They have a variational structure in 
the sense that they are obtained with a calculus of variations on functionals and there 
resolutions give explicitly the critical points of these functionals. Our aim in this paper 
is to provide them numerical schemes preserving this strong characteristic at the discrete level. 

Moreover, let us make the following important remark: since a fractional Pontryagin's system 
emerges from a fractional optimal control problem, the main unknown is then the control u. 
Consequently, the convergence of the numerical scheme constructed in Section [2] is going to 
be considered only with respect to u. 



2. Variational integrator for fractional Pontryagin's systems 

In general, fractional differential equations are very difficult to solve. One can find some 
solved examples in |25L 132), 134) using Mittag-Lefner functions, Fourier and Laplace trans- 
forms. Additionally, fractional Pontryagin's systems, as fractional Euler-Lagrange equations 
provided in pQ, present an asymmetry in the sense that left and right fractional derivatives 
are involved. It is an additional drawback in order to solve explicitly the most of fractional 
Pontryagin's systems. In this section, we then develop a numerical approach treating them. 

Nevertheless, as we have seen in Section 11.21 a fractional Pontryagin's system admits an 
intrinsic variational structure: its solutions correspond to the critical points of a functional. 
In this paper, we want to construct a numerical scheme for fractional Pontryagin's systems 
preserving at the discrete level this strong property. 

A variational integrator is a numerical scheme preserving the variational structure of a system 
at the discrete level. Precisely, let us consider a differential system coming from a variational 
principle (i.e. its solutions correspond to the critical points of a functional). Then, a varia- 
tional integrator is the numerical scheme constructed as follows: 

— firstly, one have to define a discrete version of the functional; 

— secondly one have to form a discrete variational principle on it. 

Hence, a numerical scheme is obtained and it is called variational integrator. It preserves the 
variational structure at the discrete level in the sense that its discrete solutions correspond 
to the discrete critical points of the discrete functional. Let us remind that variational 
integrators are well-developed for classical Euler-Lagrange equations in )21[ 128] and let 
us remind that we have developed a variational integrator for fractional Euler-Lagrange 
equations in [11] . In this section, we are going to construct a variational integrator for 
fractional Pontryagin's systems. 



Let us introduce the following notations available in the whole paper. Let N G N*, h = 
(b — a)/N denote the step size of discretization and T = (tk)k=o,...,N = (« + kh)k=o,...,N be 
the classical partition of the interval [a, b]. Let us assume that N is sufficiently large in order 
to satisfy the following condition: 

2h a M < 1, (cond h) 

where M is the Lipschitz coefficient of the constraint function /, see Condition ( f x Epj ). 



2.1. Reminder about discrete fractional derivatives of Griinwald-Letnikov. — 

For the sequel, we need the introduction of discrete operators approximating the fractional 
derivatives of Riemann-Liouville and Caputo. As in |13|, 115] . let us define A" and A" the 
following discrete analogous of D" and D°L respectively: 

G ^ [j^J2 a r G k-rj 

V r=0 / fc=l,...,JV 

and 

A a . ( R d)JV+l _> ( R d)7V ( 12 ) 



G ^ 7^ Yl a r G k+r 



\ r=0 / fc=o,...,JV-l 

where the elements (a r ) re M are defined by «o := 1 and 

These discrete fractional operators are approximations of the continuous ones. Indeed, 
passing to the limit h — > 0, these discrete operators correspond to the definition of the 
fractional derivatives of Griinwald-Letnikov (1867) coinciding with the Riemann-Liouville's 
ones. We refer to 1 32 1 for more details. 



Finally, according to Equation (JU), we define C A" and C A? the following discrete analogous 
of C D°_ and C -D+ respectively: 

(14) 



c^-* — 


: {R d ) N+1 - 

G H 


-> {R d ) N 


A a ■ 


(Rd)N+l _ 

G t— 


> ((AJ( G - 0K )) t ) M ^ 



and 

(15) 



Let us note that we preserve some continuous properties at the discrete level. In particular, 
G = (resp. G N = 0) implies C A°G = A^G (resp. C A^G = A%G). Additionally, in the 
classical case a = 1, these discrete fractional derivatives coincide with the usual backward 
and forward Euler's approximations of d/dt with a (—1) term in the right case: 

Vfc = 1, . . . , TV, (AlG) fc = ( c AiG) fc = Gfc ~ Gh - X (16) 



and 



Vfe = 0, . . . , N - 1, (A\G) k = ( c AlG) k = ° k ^ k+1 . (17) 
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2.2. Results concerning the discrete fractional derivatives. — In this section, we 
prove two important properties preserved from the continuous level to the discrete one. For 
the sequel, we first need the introduction of the following shift operators: 

a : (R n ) N+l — ► (R n ) N and a' 1 : (R n ) N+l — ->• (^ n ) N 

G i > ( G Wl)fe =0 ,...,JV-l G ' > i G k-l)k=l,...,N> 

(18) 
where the integer n is d or m. 

The first property is the following: considering the quadrature formula of Gauss as approxi- 
mation of the integral, we can prove the following discrete fractional integration by parts: 

Property 4 (Discrete fractional integration by parts). - Let G 1 , G 2 E (M d ) Ar+1 sat- 
isfying Gq = G 2 N = 0, then we have: 

N N-l 

hY.ic^-G'h ■ a'\G 2 ) k = hj^ a{G x ) k • ( c A«G 2 ) fc . (DFIBP) 

fe=l fc=0 

Proof -- Since G l = G 2 N = 0, we have cA^G 1 = A^G 1 and C A^G 2 = A^G 2 . Then, we 
have: 

N N-l N-lk+1 

/^(A^G 1 ), • o~\G 2 ) k = fc £(A«GVi ■ G\ = h l ~ a Y.Y. a r G Ui-r ■ G\. (19) 

k=l k=0 fc=0 r=0 

Finally, since G\ = G 2 N = 0, the following equalities hold: 

N N-l k 

h^-G'h-O-'iG 2 ), = h^Y.Y. a r G Ul-r-Gt 



k=l fc=0 r=0 

N-l N-l N-lN-r-1 



^EE^W-r^ = ^E E <*rG\ +1 -Gl +r (20) 

r=0 k=r r=0 A;=0 

iV-1 /N-k-1 \ N-l /N-k \ 

= ^- tt E G fc+l- E a rG 2 +r ) = V-o^Gl^- (Y,*rG 2 +r , 

fc=0 V r=0 / fc=0 \r=0 / 

which concludes the proof. D 

This last result is very useful for discrete calculus of variations involving discrete fractional 
derivatives, see proof of Theorem [7J Secondly let us prove the following discrete version of 
the fractional Cauchy-Lipschitz Theorem proved in [9]: 

Theorem 5 (Discrete fractional Cauchy-Lipschitz theorem) 

Let F E ^{W 1 x [a, b],R d ) satisfying the following Lipschitz type condition: 

3K eR, V{ Xl ,x 2 ,t) E {R d ) 2 x [a, b], \\F(x u t) - F(x 2 ,t)\\ < K\\xi - x 2 \\, (21) 

with h a K < 1. Then, the following discrete fractional Cauchy problem: 

C A«Q = F(Q,T) 



has an unique solution Q E (R ) 



d\N+l 



Proof. - - We are going to construct by induction the solution Q of (|22|) . Our method uses 
the classical fix point theorem concerning the contraction mappings. Indeed, let us choose 
Qo = A. Then, for any k = 1, . . . , N, Q k has to satisfy: 



fc-i 



Q k = h a F(Q k ,t k ) + Qo~Y, a r (Qk-r - Qo)- (23) 

r=l 

However, for any k = 1, . . . , N, the application h a F(-,t k ) + Qo — Ylr=i a r(Qk-r — Qo) is a 
contraction and consequently admits an unique fix point. Hence, we first construct Q\, then 
Q2, etc. By induction, we construct a solution Q of (j22|) and such a construction assures its 
uniqueness. □ 

2.3. First step of construction. — As said in introduction of this section, in order to 
complete the first step of construction of a variational integrator, we have to provide a discrete 
version of C a . In this way, let us give the following definition: 

— The elements U £ (W 11 ) +1 are called the discrete controls; 

— For any discrete control U, let Q ,a £ (IR")" +1 denote the unique solution of the 
following discrete Cauchy problem: 

c A«Q = f(Q,U,T) , 

Qo = Ae R d . l Q> 

QU,a j g ca vj ec [ the discrete state variable associated to U. Its existence and its uniqueness 
are provided by Theorem [5] and Conditions (f x lip) and (jcond h\\ ; 



— Finally, we define the following discrete cost functional: 

N 

u ^ hJ2L(Q)! ,a ,Uk,t k ). 

k=l 

Hence, we have provided a discrete version C^ to the cost functional C a . Now, the second 
step of the construction of the variational integrator is to characterize the discrete critical 
points of the discrete cost functional £^ with the help of a discrete calculus of variations. 

Let us make the following remark: such a characterization implies to be a necessary condition 
for the existence of an optimizer of the discrete cost functional C^. In fact, in this section, 
we have defined an actual discrete fractional optimal control problem. 

2.4. Second step of construction. — The second step of construction of a variational 
integrator consists in forming a discrete variational principle on £%. Precisely, we focus on 
the characterization of its discrete critical points, i.e. the elements U € (M m ) Ar+1 satisfying: 

VJ7 G (R m ) N+ \ DCl(U){U) := lim C t(U + eU) - C^(U) = Q 

With a discrete calculus of variations, we obtain the following discrete version of Lemma Q] 
giving explicitly the value of the Gateaux derivative of Cf L . 
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Lemma 6. - Let U, U € 



pm\7V+l 



. Then, the following equality holds: 



N 



DC a h {U){U) = hY J 



k=l 



-jr (Qk ' > Uk,tk) ■ Qk + -7— (Q k ' ,Uk,t k )-Uk 



dx 



dv 



(26) 



where Q E (M d ) Ar+1 is the unique solution of the following linearised discrete fractional Cauchy 
problem: 



A a _Q = %(Q U ' a ' U,T)xQ + ^(Q U ' a , U,T)xU 



(LCPg 



Oo = 0. 
Proof. - - See Appendix IA.31 



qi 



U 



This last result does not lead to a characterization of the critical points of £% yet. As in the 
continuous case, we then introduce the notion of discrete adjoint variable: for any discrete 
control U, let P U ' a £ (M. ) +1 denote the unique solution of the following shifted discrete 
Cauchy problem: 

r)J-f 

{a(Q u n,a(U),P,a(T)) 



-A1P 



P; 



N 



Or 

dx 
0. 



(^{a(Q U ' a ),a(U),a(T))'\ xP ^^ 



pu,a j g ca ri ec j i^q discrete adjoint variable associated to U . Its existence and its uniqueness 
are provided by the analogous of Theorem [5] for right discrete fractional derivative and 
by Conditions ( f x lip ) and (jcond h\\ . Let us note that, since P N ' a = 0, we can write 
A<xpU,a _ A!P U ' a . 



The presence of shift operators in the definition of the discrete adjoint variable is the 
consequence of the change of sums in the discrete fractional integration by parts (jDFIBPp 
(see Property H|). We refer to the proof of Theorem [7] for more details. We also refer to 
Remark [9] for a discussion about the presence of the shift operators. 

Finally, let us note that for any discrete control U, the couple (Q ' a , P u,a ) is solution of the 
following shifted discrete fractional Hamiltonian system: 

f)J-f 
C A«Q = —(Q,U,a-\P),T) 

^ (- HS £) 

A%P = —{a(Q),a(U),P,a(T)). 



Finally, the introduction of this last discrete element allows us to prove the following theorem: 

Theorem 7. - - Let U 6 (R m ) Ar+1 . Then, U is a discrete critical point of £% if and only if 
(Q u,a , U, P U ' a ) is solution of the following shifted discrete fractional stationary equation: 



^(Q,U,a-\P),T) 



0. 



(aSE? 



11 



Proof. -- Let U,U £ (W n ) N+1 . Prom Lemma El we have: 



A' 



h- 1 DCt{U){U) = Y J 



fe=i 



(Q"' a ,U k ,t k )+(^(QV' a ,U k ,t k )\ xo- l (P u ' a )k 



■Qk 



-^((%(Qk^^k,t k )) T xa-^P^^,) ■Q k + J2^(Q k J > a ,U k ,t k ).U k . (27) 

fc=i V ^ ' J fc=i 

Then: 

h^DCUUW) = J2 (c& a + P u ' a )k ■ <?{Q) k - J2 /(Q? ,a > ^> **) x ^ ) ' ff " 1 ( pU,a )* 

fc=o fc=i ^ x ' 

N BT 
+ Y,o-(Qk' a ,Uk,t k )-U k . (28) 

From the discrete fractional integration by parts (jDFIBPp (see Property 2]) , we obtain: 

h- l DC a h {U){U) = J2 ((c^Q)k - j-(Q%' a , Uk,t k ) x Q k ) ■ a~\P u ' a )k 
k=i ^ X ' 

^ ar 

+ X)^(Q? ,a ,^,**)-^- (29) 



fc=i 



Since Q is solution of (TLTTPlO), we have 



Q 



N rat \ N bt 

h- l DC a h {u)(V) = Y, [j-(Q%' a ,Uk,t k )xu k )-°~\P u nk + Y.-^ (Q^ a ,u k ,t k )-u k 

k=\ ^ V ' fe=l 



fc=l 
Finally: 



r as 



D£%(U)(U) = hJ2-Q- {Q l i ,a ^ k ^-\P u ' a )k,t k )-Uk. (30) 

fc=l 

The proof is completed. D 

Finally, from Theorem [71 we obtain the following result leading to the variational integrator 
constructed: 

Corollary 8. -- £% has a discrete critical point if and only if there exists (Q, U, P) £ 
(M. d ) N+1 x (M"*)^" 1 " 1 x (IR^)^" 1 " 1 solution of the following shifted discrete fractional Pontryagin's 
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system: 



c A a _Q = —{Q,U,a~ 1 (P),T) 
A%P = ^-(a(Q),a(U),P,a(T)) 



¥L{q,u,*-Hp)>t) 



(aPS^ 







I (Qo,Pn) = (A,0). 
In this case, U is a discrete critical point of C^ and we have (Q,P) = (Q ' a ,P ' a ). 

Let us note that (jcrPS^ ) is made up of the shifted discrete Hamiltonian system (crHS^), the 



shifted stationary equation (jcSE^h and initial and final conditions. 



Hence, we have constructed the variational integrator (|crPS?p for the fractional Pontryagin's 



system (|PS Q p . It is then a numerical scheme for (|PS"P preserving its variational structure in 
the sense that the discrete solutions U obtained correspond to the discrete critical points of 
the discrete version C^ of C a . 



Remark 9. - Let us note that the variational integrator (cPS^) does not correspond with a 
direct discretization of (|PS a p as it is done in |13j . There is an emergence of shift operators 
caused by the conservation at the discrete level of the variational structure. However, it is 
proved that the use of shifted numerical schemes allows to obtain more stability for some 
fractional differential equations, see |29L 130] . 

Remark 10. - Let us remind the following remark: since a fractional Pontryagin's system 
emerges from a fractional optimal control problem, the main unknown is then the control u. 
Consequently, the convergence of the variational integrator (aPS^;) is going to be considered 



only with respect to u. Let us note that the value of Uq does not take place in the variational 
integrator (cPS^): it is a free value. Nevertheless, this is totally coherent with the fact that 
this value does not take place neither in the definition of C^ . Hence, in the following examples 
in Section 0, the error between an exact solution u of (|PS a p and a numerical solution U 
obtained with (ctPS?) is going to be evaluated on ||u(ifc) — Uk\\ f or k G {1, . . . , N} only. 



2.5. Link with the discrete fractional Euler-Lagrange equation 

the constraint function f(x,v,t) 



- Let us take 
In this case, applying Corollary 



v satisfying <\f x lipP 
[HI we know that there exists a critical point of £% if and only if there exists a solution 

of the shifted discrete fractional Pontryagin's 



{Q,u,P) e (irr +1 x (R m ) 

system (krPS^p here given by: 



m\N+l 



N+l 



A-Q = u 



A«P 
dL 



Ov 
I (Qo,Pn 



^(a(Q),a(U),a(T)) 

(Q,U,T) + a~ 1 (P) = 
(A,0). 



(31) 
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In the affirmative case, it implies that Q is a discrete solution of the following discrete frac- 
tional Euler- Lagrange equation: 



g(Q, c A«Q,r) + A ? (^ 



— (Q, C A^Q, T) + A%[ — (Q, C A! Q, T) ) = 0. (EL£ 



Finally, according to our works in [11] . we then obtain that Q is a critical point of the 
following discrete fractional Lagrangian functional: 



N 



Q^/ l ^L(Q fe ,( c A^Q) fc ,i fc ). (32) 



fc=i 



We refer to |llj for more details concerning discrete fractional Euler-Lagrange equations. 



3. Numerical tests 

In the following numerical tests, according to Remark [10\ we are going to give graphic rep- 
resentations only of discrete solutions U and the study of the convergence of the variational 
integrator (uPS^) is only going to be evaluated on the convergence of the discrete control to 



the continuous one. 



3.1. The linear-quadratic example. — Linear-quadratic examples are often studied in 
the literature because they are used for tracking problems. The aim of these problems is 
to determine a control allowing to approach as much as possible reference trajectories, |36[ 
Part 1.4, p. 49]. In this section, we study such an example, [16} Part 4.4.3, example 3, p. 53]. 
More generally, a quadratic Lagrangian is often natural (for example in order to minimize 
distances) and even if the constraint functions are frequently non linear, we are often leaded 
to study linearised versions. 

Let us choose d = m = A = 1 and [a, b] = [0, 1]. Then, let us take the following quadratic 
Lagrangian and linear constraint function: 



L: R 2 x [0, 1] — ► R and / : K 2 x [0, 1] — ► R (33) 

(x,v,t) i — > (x 2 + v 2 )/2 (x,v,t) i — ► x + v. 
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Let us give the graphic representations of the numerical solutions U given by (crPS^) for 
N = 500 and for a = 1, 3/4, 1/2, 1/4: 



Discrete solution U 




We have seen in [S] that the fractional Pontryagin's system (|PS"P is explicitly solved only in 
the classical case a = 1 and we obtained the following unique critical point of C : 



VtG [0,1], n(t) 



cos 



h(V2) 



i? 



■ sin 



h(v^t) 



sin 



h(V2) 



i? 



cos 



h(>/2t), 



(34) 



where R = y/2 cosh(\/2) — sinh(\/2). Hence, we can only test the convergence of the variational 
integrator (jaPS?) for a = 1. We give the following graphic representing the logarithm of the 



error max (\u(tk) — Uk\, k = 1, . . . , N) versus the logarithm of h and the identity function for 
comparison: 

alpha=1 



-2.5 

-3 
-3.5 

-4 
-4.5 

-5 
-5.5 

-6 
-6.5 

-7 




-6.5 



-5.5 -5 

log(h) 



-4.5 



-3.5 



In this example with a = 1, the convergence seems then obtained with order 1. Neverthe- 
less, we do not know the exact solution of (|PS a p in the strict fractional case < a < 1. 
Consequently, we can not study the behaviour of the error in this case. 
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3.2. A solved fractional example. — In this section, we are going to compute (crPS^) 
in the framework of an example solved in the strict fractional case in the sense that we know 
explicitly the unique critical point u of £% for any < a < 1, see [9]. Consequently, for this 



example, we can test the convergence of the variational integrator (<rPS?) for any < a < 1. 



Then, let us choose d = m = A = 1 and [a, b] 
Lagrangian and linear constraint function: 



[0,1]. Then, let us take the following 



L: R 2 x [0, 1] 



R and / : M 2 x [0, 1] 

(l-t)x + (v 2 /2) (x,v,t) 



(35) 



X + V. 



Let us give the graphic representations of the numerical solutions U given by (<rPS^) for 
N = 500 and for a = 1, 3/4, 1/2, 1/4: 



Discrete solution U 




-0.5 

-1 
-1.5 

-2 
-2.5 

-3 
-3.5 

-4 
-4.5 



/ 


- 




1 U 4 















0.2 



0.6 



As we have seen in [9], the fractional Pontryagin's system (|PS"j) is explicitly solved for any 
< a < 1 and we obtained the following unique critical point of C a : 



V£E [0,1], u(t) 



tr +i E a , a+2 ((i-tr), 



(36) 



where E Qia+ 2 is the Mittag-Lemer function with parameter (a, a + 2). Let us the convergence 



of the variational integrator (erPS^) for any < a < 1. We give the following graphics 



representing the logarithm of the error max (|u(£fc) — Uk\, k = 1, . . . , N) versus the logarithm 
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of h and the identity function for comparison for a = 1, 3/4, 1/2, 1/4: 



alpha=1 




alpha=3/4 




-5.5 -5 

log(h) 



alpha=1/2 



alpha=1/4 




-5.5 -5 

log(h) 




-5.5 -5 

log(h) 



For this example, the convergence seems then obtained for any a = 1, 3/4, 1/2, 1/4 and still 
with order 1. Hence, the graphics obtained in these Sections 13.11 and 13.21 make us confident 



with respect to the quality of (crPS?) both in the classical and strict fractional cases. 



4. A discrete fractional Noether's theorem 



Fractional Pontryagin's systems (|PS Q jl are very difficult to solve explicitly, see example 
in Section 13. 11 Consequently, a deviously way in order to get informations on the exact 
solutions is to derive a constant of motion, i.e. functions which are constant on each solution. 
Such conservation laws allow to obtain many informations in the phase space for example 
or to integrate the equation by quadrature. In [9], we prove a fractional Noether's theorem 
giving the existence of an explicit conservation law for fractional Pontryagin's systems (|PS"D 
exhibiting a symmetry. Let us remind that this result is based on a preliminary result proved 
by Torres and Frederico in |17|, 118] . 

In this section, we study the existence of discrete conservation laws for shifted discrete 
fractional Pontryagin's systems (<tPS?). Precisely, following the same strategy, we introduce 



the notion of discrete symmetry for such systems and prove a discrete fractional Noether's 
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theorem providing an explicit computable discrete constant of motion. Let us note that this 
work is strongly inspired from our study in |10j where we have provided a discrete fractional 
Noether's theorem for discrete fractional Euler-Lagrange equations admitting a discrete 
symmetry. 

We first review the definition of a one parameter group of diffeomorphisms: 

Definition 11. -- Let n G N*. For any real s, let <p(s,-) : R n — > W 1 be a diffeo- 
morphism. Then, $ = {(f>(s, -)} s eU is a one parameter group of diffeomorphisms ofW 1 if it 
satisfies: 

1. 0(O,-)=Id R »; 

2. Ms, s' £ R, (j)(s, •) o <j)(s', ■) = (p(s + s', •); 

3. 4> is of class < ^ 72 . 

Usual examples of one parameter groups of diffeomorphisms are given by translations and 
rotations. The action of three one parameter groups of diffeomorphisms on an Hamiltonian 
allows to define the notion of a discrete symmetry for a shifted discrete fractional Pontryagin's 



system (<rPS?): 



Definition 12. -- Let <3?j = {4>i(s, -)} se iR ; for i = 1,2,3, be three one parameter groups of 
diffeomorphisms of R d , M. m and M. d respectively. Let L be a Lagrangian, f be a constraint 
function and H be the associated Hamiltonian. H is said to be c A" -invariant under the 
action of (3>i)i=i,2,3 if it satisfies: for any (Q, U, P) solution of (crPS^p and any s£R 



iJ^i( S ,Q),0 2 ( S ,J7),0 3 ( S) a- 1 (P)),Tj -(f> 3 (s,a-\P)) ■ c A°L(<l>i{s,Q) 

= H(Q, U, <r-\P),T) - a~\P) ■ c A a _Q. (37) 

From this notion, we prove the following Lemma: 

Lemma 13. Let L be a Lagrangian, f be a constraint function and H be the associated 
Hamiltonian. Let us assume that H is C A™ -invariant under the action of three one parameter 
groups of diffeomorphisms (3>i)i=i2,3- Then, the following equality holds for any solution 



(Q,U,P) solution of flgPSgp : 

^(0, Q) ■ a~\A%P) - C A° (^(0, Q)\ ■ a~\P) = 0. (38) 



Proof. - - Let us differentiate (|37|) with respect to s and let us invert the operator C A" and 
d/ds. Taking s = 0, we finally obtain: 

dH . . d&i , _. dH . . d(t>2 , TTN OH, . 863, i,„ N \ 

- ^(0, ^(P)) ■ cA a _Q - a-\P) ■ C A^ (^(0, Q)\ = 0, (39) 
where *= (Q, U, a' 1 (P),T). Since (Q,U,P) is solution of (|ffPSgp, we obtain (15511 . □ 



Let us note that this last result corresponds to the discrete version of the result proved by 
Torres and Frederico in |17|, I18j . Let us remind that our aim is to provide an explicit discrete 
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constant of motion for shifted discrete fractional Pontryagin's systems (erPS^) exhibiting a 
discrete symmetry. Our result is based on Lemma [T3l and on the following implication: 

vggm jv+1 , aLg 







3c G R, V/c = 0, . . . , N, G k = c. 



(40) 



Namely, if the discrete derivative of G vanishes, then G is constant. Consequently, our aim 
is to write the left term of (|38p as an explicit discrete derivative (i.e. as Al of an explicit 
quantity) . In this way, we are going to use a discrete transfer formula as it is done in |10] 
for discrete fractional Euler-Lagrange equations admitting a discrete symmetry. 

Nevertheless, we have first to introduce some square matrices of length (A + 1). First, 
B\ := Idjv+i and then, for any r E {2, . . . , A}, the square matrices B r G Mn+i defined by: 

Vi, j =0,...,N, (B r )ij : = £{i<j<]V-l}£{l<j<iV-r}£{0<i_j<r-l} - <5{ J=0 }<5{ r <j}, (41) 

where 5 is the Kronecker symbol. Secondly, we define the square matrices C r G Mn+i by: 

Vr = l,...,A, Vi,j = 0,...,A, (C r )ij := <5 {r < i} <5 0=0} . (42) 

Finally, we define the square matrices A r E 7Wat+i by: 

Vr = 1, . . . , TV, A r := a r 5 r + 0?C r , (43) 

where /?" = XH-=o a fc- Examples of matrices A r G A^at+i for N = 5 are given in Appendix 
[A3 

Lei G 1 , G 2 G (R^)^ 1 satisfying G 2 N = 0. 



Lemma 14 (Discrete transfer formula). 

Then, the following equality holds: 



N 

G 1 ■ a-\A1G 2 ) - (cA^G 1 ) • ^(G 2 ) = /i^A^i, x (G 1 • ^(G 2 )) 



r=l 



Proof — See Appendix IA.4I 

Consequently, combining Lemmas [13] and [131 we prove: 



(44) 
□ 



Theorem 15 (Discrete fractional Noether's theorem). - Let L be a Lagrangian, f be 
a constraint function and H be the associated Hamiltonian. Let us assume that H is C A°- 
invariant under the action of three one parameter groups of diffeomorphisms ($i)i=i,2 3- Then, 



the following equality holds for any solution (Q,U,P) of (uPS^): 

N 



5>x(^(0,Q).^(P)) 

r=l ^ ' 



0. 



(45) 



According to Equation (|40p . this theorem provides a discrete constant of motion for any shifted 



discrete fractional Pontryagin's systems (uPS^) exhibiting a discrete symmetry. Moreover 



this discrete conservation law is not only explicit but also computable in a finite number of 
steps. Let us see a concrete example: 



Example 16. Let us consider d = m 
following linear constraint function: 



L 



i 2 xR 2 



x [a, b] 
(x,v,t) 



2, the following quadratic Lagrangian and the 
and f : 



^ 2 xR 2 



klr + lbll)/2 



x [a, b] 
(x,v,t) 



X + V. 

(46) 
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Then, we consider the three one parameter groups of diffeomorphisms given by the following 

rotations: 

fa: IxR 2 — ► R 2 (47) 

cos(s6>i) — sin(s6>j) \ / x\ 
sin(s^j) cos(s9j) J \ x 2 

for i = 1, 2, 3 and where Q\, 9 2 E R and #3 = —9\. With these parameters, one can prove that 
the Hamiltonian H associated to L and f is c A" -invariant under the action of ($1)1=1,2,3 • 
Consequently, the fractional Pontryagin's system (crPS^) admits a symmetry and then admits 



(s,x 1 ,x 2 ) 



an explicit discrete conservation law given by the discrete fractional Noether's Theorem \15[ 



We choose A = (1, 2), N = 100 and 6 1 = 6 2 = -6 3 = 1. Let us compute ( jcrPSgp for a = I, 
3/4, l/2 ; 1/4. Then, we denote Q = (Q l ,Q 2 ) and P = {P l ,P 2 ) the discrete solutions 
obtained and we denote G = d(j)i/ds(0,Q) = (—Q 2 ^ 1 ). We are then interested in the value 
of: 



N 

E 

r=l 



A r x (G-<r r - 1 (P)). 



(48) 



Let us see the graphics obtained by the computation of (aPS^) and by the computation of the 
quantity given in Equation (f48]) for a = 1, 3/4, l/2 7 1/4; 



Discrete solutions (Q1 ,Q2) and (P1 ,P2), alpha=1 



Discrete conservation law on the solutions (Q1 ,Q2) and (P1 ,P2), alpha=1 





Discrete solutions (Q1 ,Q2) and (P1 ,P2), alpha=3/4 




Discrete conservation law on the solutions (Q1 ,Q2) and (P1 ,P2), alpha=3/4 



- constant of motion 
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Discrete solutions (Q1 ,Q2) and (PI ,P2), alpha=1/2 



Discrete conservation law on the solutions (Q1 ,Q2) and (PI ,P2), alpha=1/2 










| constant of motion 


0.8 




0.6 




0.4 


- 


0.2 


- 









0.2 




0.4 


- 


0.6 


- 


0.8 





Discrete solutions (Q1 ,Q2) and (P1 ,P2), alpha=1/4 




Discrete conservation law on the solutions (Q1 ,Q2) and (P1 ,P2), alpha=1/4 



- constant of motion 



As expected from Theorem \ 151 we obtain discrete constants of motion for this discrete frac- 
tional Pontryagin's system (oFS?) admitting a discrete symmetry and for any a = 1, 3/4, 



1/2, 1/4. In this specific example, the constant obtained is zero. 



A 

Appendix A 

A.l. Result of stability of order 1. — In this section, we use the notations and definitions 
given in Sections 12.31 and 12.41 Let us prove the following Lemma: 

Lemma 17. Let U, U € (R ?n ) iV+1 . Then, there exists a constant C\ > such that: 

V|e|<l, Vk = 0,...,N, \\Q^ +£ °' a -Q^ a \\ <C x \e\. (49) 

Proof. — First, let us prove by induction that: 

(50) 



Vfc = 0, . . . , N, 3R k > 0, V|e| < 1, \\Q^ +eU ' a 



W k || S JrCk\£\- 



For k = 0, let us take Rq = since Q 



Qq'° = A for any |e| < 1. Let k € {1, ... , N} 



and let us assume that the result (|5Q|) is satisfied for any n = 0, . . . , k — 1. Let us prove that 
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the result ([50|) is then satisfied for n = k. Since Q u,a (resp. Q u+eU,a ) is solution of (CPg ) 
associated to U (resp. to U + et/), we have: 



Q U k ' a = h a f(Q"' a , U k ,t k ) +A-J2 <*r(Q%l° ~ A) 



(51) 



r=l 



and for any |e| < 1: 



Q U+eU,a = haf{Q U+eU, a 



Consequently, for any |e| < 1: 



,U k + eU k ,t k )+A-^2otr(Q 



U+eU,a A , 
k-r A >- 



(52) 



r=l 



\\Q k +e "' a - Q k ' a \\ < h a \\f(Q" +e °' a ,U k +eU k ,t k ) - f(Q^ a ,U k) t k )\\ 



Y,\ a r\\\Q U k X £tJ ' a -Q U k l a rl (53) 



r=l 



and then, with the induction hypothesis, we have for any \e\ < 1: 

\ lQ U+eU,a _ Q U,a^ < fc«||/(Q^.«, ^ + e ^, tfc ) - f(Q%' a ,U k + eU k ,t k )\\ 

k 

+ h a \\f(Q^ a ,U k + eU k ,t k ) - f(Q^ a ,U k ,t k )\\ +Y,\a r \Rk-r\e\. (54) 



r=l 



Finally, using Condition (/ x lip) and a Taylor's expansion of order 1 with explicit remainder, 
we prove: 



V|e|<l, \\Q k ~Q k \\<h a M\\Q k - Q k ' \ 



+ h a \e\ 



%{Q u k ' a ,e,t k )*u k 



+ y^l a r-|-Rfc-r|g| i (55) 



r=l 



where £ e G [£4,£4 + e^fc] C [— Mi,Mi] m with Mi > independent of |e| < 1. Hence, since 
df/dv is continuous, we can conclude that there exists M 2 > such that: 



vkl < l, 



9f , n U,a 

dv {Qk 



^-(QY' a ,e,t k )xu k 



<Mo 



(56) 



Consequently, since h a M < 1 from Condition (jcond h\i . we have: 



V|e|<l, \\Q? +e0 ' a -QY' a \\< 



c 



l-h a M 



h a M 2 +J2Rk-r\ar\) ■ 



(57) 



r=l 



We then define R k := (h a M 2 + J2r=i R k-r\a r \) / (I - h a M) independent |e| < 1 which con- 
cludes the induction. To complete the proof, we have just to define C\ = m&x{R k , k = 
0,...,N}. D 
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A. 2. Result of stability of order 2. — In this section, we use the notations and definitions 
given in Sections 12.31 and 12.41 We prove the following Lemma with the help of Lemma [TTl 



Lemma 18. -- Let U, U G (l™)^ 1 . Then, there exists a constant C > such that: 

V|e| <1, Vk = 0,...,N, \\QV +£ °> a -QV> a -eQ k \\<Ce 2 , (58) 

where Q is the unique solution of the following linearised discrete fractional Cauchy problem: 



C A« Q = ^(Q U ' a , U,T)xQ + ^(Q U ' a , U,T)xU 



(LCP 



Its existence and its uniqueness are provided by Theorem [3| and Conditions (f x lip) and 
(jcond h\f . 



Proof. - - We proceed in the same manner that for Lemma [T71 Let us prove by induction 
that: 



\/k = 0,...,N, 3R k >0, V|e|<l, \\Q k +eU ' a - Q^' a - eQ k \\ < R k e 2 . 



(59) 



For k = 0, let us take Rq = since Q ,a = Q ,a = A for any |e| < 1 and Qo = 0. Let 
k G {1, . . . , N} and let us assume that the result (|5U|) is satisfied for any n = 0, . . . , k — 1. 
Let us prove that the result ([59]) is then satisfied for n = k. Since Q Ua (resp. Q u + eU ^ \ s 



solution of (LCP^) associated to U (resp. to U + ell) and since Q is solution of (LCP^), 
we have with a Taylor's expansion of order 2 with explicit remainder: 



= 1, Qk 

1 



dx 



U+eU,a _ Q U, a _ £Qk = h a^_ {Q U, aj jj^ x ( Q U+eU,a _ Q U,a _ ^ 



+■ h a { £ w(& & i fc )(Qr £t/,a - C a > ^ °) 2 - E «^ + 



' k—r 



^5j._j. £{°ik—r)i 



r=l 



(60) 



where: 

- a e [Q^ a ,Q^' a ] C [||Q^' a || - Ci,||Q?' a || + Ci] from Lemma El Then, £f G 
(-Mi,Mi] d with Mi > independent of |e| < 1; 

- £f G p7 fc , C/ifc + eU k ] C [-M 2 , M 2 ] m with M 2 > independent of |e| < 1. 

Since V 2 /(-, -,t k ) is continuous, we conclude that there exists M3 > such that: 



Vld < 1, 



v 2 /(ef, ei, t fc )(Qr £C,Q - ^' a > ^> 0) 2 



S -M 3 (J|C^ fc - Q k || +2||y fc _ Vfc llll e ^A;|| + \\ £U k\\ j 



Hence, from Lemma [TTl there exists M4 > such that: 
Vlel < 1, 



V 2 /(£f , CI, t k )(Q" +eU ' a - Q^ a , eU kl 0) 2 



< M 4 e 2 



(61) 



(62) 
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From Equality (j60|) and Condition (f x lip), Inequality ()62p and the induction hypothesis, we 
obtain: 

V | £ | < 1, ||Q™' Q - Q^ a - eQ k \\ < 2h a M\\QV +£ °' a - Q^ a - eQ k \\ 



+ ^T\a r \R k - r e 2 + h a M 4 e 2 . (63) 



r=l 



Finally, since 2h a M < 1 from Condition (|cond h\\ , we have: 
V|e|<l, \\Q U k +£0 > a -Q U k > a -eQ k \\< 



1 - 2h a M 



h a Mi + J2R 



k — r \Q r \ 



(64) 



r=l 



We then define R k := {h a M± + Yl r =i ^fc-r|ov|)/(l — 2h a M) which concludes the induction. 
In order to complete the proof, we just have to define C := max{R k , k = 0, . . . , N}. □ 

A. 3. Proof of Lemma O — In this section, we prove Lemma and consequently, we use 
notations and definitions given in Sections 12.31 and 12. 41 



Let U,U G (M m ) 7V+1 and Q e (IR^) 7 ^ 1 the unique solution of (ILCP^p. From Lemma [HI 
we have: 



Vk = 0,...,N, V|e| < 1, Q^ +£ ^ Q = QJ^' a + eQ fc + H £ k , 



(65) 

where for any k = 0,...,N and for any |e| < 1, \\H k \\ < Ce 2 . In particular, there exists 
Mi > such that: 



•.U.a ^U+eU.ai 



Vk = 0,...,N, V|e|<l, [Q^°,Qp u "*]c -M^Mi] 



In the same way, there exists M 2 > such that: 



VA; = 0, . . . , N, V|e| < 1, [Z7*, C/ fc + eL7 fe ] C [-M 2 , M 2 



(66) 

(67) 



We have: 



N 

V|e| < 1, ££(17 + efr) - ££(J7) = h £ ^(Q^^, 17* + ^*, **) - L(Q^' a , 17*, **)] ■ (68) 



*=i 



With a Taylor's expansion of order 2 with explicit remainder, we have for any |e| < 1 and 
any k = 0, . . . , N: 



'^^+^,i fc )-L(Q^^^,t fc )-e^(Q^^^,t fe ).Q fe -^(Q^ Q ,C/ fc ,i fc )-L7, 



)U,a 



■JJ,a 



< 



dL [Q U k >\U k ,t k ).Hi\ + \^ 2 

■JJ,a ^U+eU,a 



^(Q^ a ,U k ,t k )-H e k + l -V 2 L{£l&,t k ){eQ k + H e k ,eU k ,ti)\ (69) 



where £f G [Q k ' a ,Q k +£U ' a ] C [-Mi,Mi] d and £§ € [*7*, £/* + rf/*] C [-M 2 ,M 2 ] m . Since L 
is of class < ^ 72 , we obtain easily that there exists M3 > such that for any |e| < 1 and any 
fc = 0,...,JV: 

L(QV +£ °> a , U k + eU k , t k ) - L(QV' a , U k , t k ) 

BI BT 

-e-^(Qk' a ,U k ,t k )-Q k -e—(Q^' a ,U k ,t k )-U k <M 3 e 2 . (70) 
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Consequently, we have for any < |e| < 1 and any k = 0, . . . , N: 
L(QV +£ °' a , U k + eU k , t k ) - L(Q^ a , U k ,t k ) 



—-ft-{Qk' ' Ukitk) ■ Qk — Tj— (Q k ' ,u k ,t k )-u k 



<M 3 e. (71) 



Hence: 



N 



lim Ci(U + eU)-C t (U) =hE 



e^O 



fc=l 



dL 



-pr- (Q k ' ,U k ,t k ) ■ Q k + —(Q k ' ,U k ,t k )-U k 



(72) 



The proof is completed. 



A. 4. Proof of Lemma 1141 — In this section, we use the notations and definitions given 
in Section SJ Let us prove Lemma [Til 



Let G 1 , G 2 e (M. d ) N+1 satisfying G% = 0. First, let us denote for any k = 1, . . . ,N: 



X k := h a G\ ■ a-HKG 2 )k - (cA1G% ■ o~\G 2 ) h 



(73) 



Our aim is to write X as an explicit discrete derivative {i.e. as Al of an explicit quantity). 
We have for any k = 1, . . . , N: 



'N+l-k 



X k = G\-\ Y, «rG2 +r ,_ 1 )-(^a r (Gt r -Gi)j-Gt 1 =a 1 / l (A^G 1 .G 2 ) fc + n + Z fc , 

(74) 



r=0 



^r=0 



where for any k = 1, . . . , N: 



Yk '-~ I 7 , °r j ^0 " ^fe-1 - Pk ^0 " ^fc-1 



and 



2fc :- 



'N-k 



@k " I 7 , "r^fc+r-l - / v a rG k _ r ■ G k _ 1 



,r=2 



vr=2 



(75) 



(76) 



Our aim is then to write Y" and Z as explicit discrete derivatives. We then define for any 
i = 0,...,N, Vi := hY^ r= \Y r and W{ := /&X)j - =i^r Hence, we have A^V = Y and 
A\W = Z and then, X = Al(ai/iG 1 • G 2 + V + W). Our aim is then to explicit V and 
W. For any i = 0, . . . , N, we have: 



JV iv 



y, = /i J>?G$ • G 2 _ x = fcJ^Gj • a r - 1 (G 2 )o = hj^ Y,ffC r (i,j)G] ■ a r -\G 2 ) r (77) 

r=l j=0 



r=l 



r=l 
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For any i = 0, . . . , N, we have: 

j=l L \r=2 / \r=2 / 

j JV-j i j 

j=l r =2 j=2 r=2 

j AT— j i i 

j=l r=2 r=2 j=r 

i N—j i i—r 

j=l r=2 r=2 j=0 

The following equality holds for any r = 2, . . . , N and any i, j = 0, . . . , N: 



(78) 
(79) 
(80) 
(81) 



^{l<J<i}^{2<r<7V-.j} ~~ d{2<r<i}<){0<j<i-r} 

= ${l<i<N-l}${l<j<N-r}3{0<i-j<r-l} ~ <%=()} <5{ r <i}- (82) 

Consequently we have for any % = 0, . . . , N: 

N N 

Wi = hJ^Yl a r Rr (*' fi G ) • fjr_1 ( G2 )j • ( 83 ) 

r=2 j=0 

Finally, we have for any i = 0, . . . , N: 

N N 

ai hGJ ■ G\ + Vt + Wi = hJ2J2Mi,J)Gj ■ ° r ~HG 2 )r (84) 



r=l j=0 



Finally, the following equality holds: 



N 

X = hA 1 _\j2A r x (G l ■ <r r -\&)) 



r=l 



(85) 



which concludes the proof. 



Now, let us see some examples of matrices A r £ .Mtv+i for N = 5: 

/ai00000\ / 00000\ 



A, 



(8f ai 

/3f ai 

/3f ai 

/3f ai 

\ /3f ai / 

/ o ooooo\ 

a 3 

a 3 a 3 

/?3 — a 3 as as 

P^ -as a 3 

\/3f-a 3 0/ 



-4, 



a 2 

$2 ~~ a 2 a2 ° 2 ^ ^ ^ 

@2 ~ a 2 02 «2 

/3f - a 2 a 2 

y/3n?-a 2 00/ 

/ 00000\ 

a 4 

a 4 

a 4 

/3f - a 4 a 4 

\/3£-a 4 0/ 
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and 



/ 00000\ 









\ Pg-a 5 0/ 
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